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ABSTRACT 

This paper reports the detection and the measurements of occultations of 
the two transiting hot giant exoplanets Kepler-5b and Kepler-6b by their parent 
stars. The observations are obtained in the near infrared with Warm-Spitzer 
Space Telescope and at optical wavelengths by combining more than a year of 
if e|}/er photometry. The investigation consists of constraining the eccentricities of 
these systems and of obtaining broad band emergent spectra for individual plan- 
ets. For both targets, the occultations are detected at 3 a level at each wavelength 
with mid-occultation times consistent with circular orbits. The brightness tem- 
peratures of these planets are deduced from the infrared observations and reach 
Tspitzer=1930 ± 100 K and Tspiteer=1660 ± 120 K for Kepler-5b and Kepler-6b 
respectively. We measure optical geometric albedos Ag in the Kepler bandpass 
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and find Ag = 0.12 ± 0.04 for Kepler-5b and Ag = 0.11 ± 0.04 for Kepler-6b 
leading to upper an limit for the Bond albedo Ab < 0.17 in both cases. The 
observations for both planets are best described by models for which most of the 
incident energy is redistributed on the dayside, with only less than 10% of the 
absorbed stellar fiux redistributed to the night side of these planets. The data for 
Kepler-5b favor a model without a temperature inversion, whereas for Kepler-6b 
they do not allow distinguishing between models with and without inversion. 

Subject headings: planetary systems — eclipses — techniques: photometry 



1. Introduction 



The highly irradiated transiting hot-Jupiters currently provide the best opportunities 
for studying exoplanetary atmospheres in emission, during planetary occultations, when the 
exoplanets pass behind their parent stars. Emission from transiting hot-Jupite rs (HJs) have 
been studied theoretically prior to any detection of light from an alien world ( ISeager et al. 
2000l : ISudarsky et al.ll2000l ). The first exoplanetary emi s sions were first detec ted from space 



in the infrared wavelengths (ICharbonneau et al. 



cently in the optica l bandpass (lAlonso et al 



2005 



2009a 



Deming et al 



Snellen et al 



2009 



Alonso et al.l l2010l ) using the CoRoT space telescope. Notably, iRowe et al.l ( 12006 



2005) and more re- 



Borucki et a 



2009 



20081) 



use the Microvariablity and Oscillations of Stars {MOST) telescope to place a very strin- 
gent upper limit on the optical depth of the occultation of HD 209458b. Recently, ground 
based spectrophotometry has also permitted observations of HJs in the near infrared (e.g . 
Sing k L6pez-Moralesl[2009l : llopez-Morales et al.l[201o[ Icillon et al]l2009[ ICroU et al.lboioh . 
The planet-to-star fiux ratio is in the order of ten times higher at infrared wavelengths, due 
to the thermal emission from the pl anet, compared to optical wavelengths, domain where the 
starlight is refiected by the planet (ISudarsky et al.ll2003l : iFortney et al.ll2005t iBurrows et al. 



20051 : ISeager et al.l l2005l : iBarman et al. 



20051 ). Consequently, only a few objects can be de- 



tected in the visible, those for which we will have enough observations to combine during 
occultations to improve the signal-to-noise and reveal the eclipse event. However, these ob- 
jects are very valuable for understanding the energy budget of HJs and, to some extent, 
for comparative exoplanetology. The HJs detected by CoRoT and Kepler ar e particularly 
interesting since their oc cultations can be observed in the optical bandpass (ISnellen et al. 
2OO9I : iBorucki et al.ll2009l ). 



Obtaining multiple wavelength measurements of the relative depths of planetary occul- 
tations is fundamental to understanding the energy budget of these objects ( ISudarsky et al. 
2OO3I : iBurrows et al.ll2005[ ). The thermal structure of HJ is likely to be infiuenced by strong 
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stellar irradiation. For example, high altitude optical absorbe r s can lead to atmospheric 



temperature inversions (IHubeny et al.ll2003l : iFortnev et al- 



as observed to a wide range of irradiation levels (e.g. [Harrington et al. 



2008; 



Burrows et al. 2007.12008) 



2007 



f 



Knutson et al. 



20081 : iMachalek et all l2008l : iTodorov et al.ll2010f ). Although the nature of the absorber is 



subject to debate, titanium oxide (TiO) and vanadium oxide (VO) are amongst the best 
molecular candidates, since they are expected to be present in these hot atmospheres an d 
that they are strong absorbents of the incident visible stellar flux (IFortney et al.l l2008l ). 
Furthermore, non-equilibrium chemistry, due to photochemistry or vertical mixing, could 
potentially lead to strong UV or optical absorbers being found in the upper atmosphere of 



these planets. This can also play a role in inverting the thermal structure (IBurrows et al. 



20081 : IZahnle et al.ll2009l ). although degeneracies betwe en thermal inversions and che mical 



composition are known to exist in atmospheric models (IMadhusudhan &: Seagerll2009l ). In- 
terestingly, the presence of high altitude absorbers in HJ at mospheres can potentially be 



revealed with transmission spectroscopy (jPesert et al. 



2008h. Emission spectroscopy ob- 



serva ti ons help also t o prob e the atmosphere of HJs (IKnutson et al.l l2008l : iMachalek et al. 



20091 ) . ISpiegel et al.l (l2010al ) and ISpiegel fc Burrowsl (l2010b[ ) demonstrate how the Kepler- 



hand optical flux from a hot exoplanet depends on the strength of a possible extra optical 
absorber in the upper atmosphere. Therefore, multi-wavelength observations are necessary 
to determine the energy budget of these exoplanets. Measurements of the optical and near 
infrared emergent fluxes from the planet allow us to derive the Bond albedo of these objects, 
to learn about the energy distribution and to address whether or not an atmospheric thermal 
inversion can be revealed. In the case of HD 209458b, the MQiS'T observ ations show that the 
planet has a very low albedos at optical wavelengths (IRowe et al.ll2008l ). which rules out the 
presence of br i ght re flective clouds in this exoplanet's atmosphere. In a similar framework. 



Deming et al.l (120101 ) combine occultations from Co Ro T and from various spa ce and ground- 



based infrared measurements of CoRoT-1 and 2 ( Gillon et al. 20091 . boiol ). to probe the 
atmospheric structures of these HJs. The multi-wavelength observations show that CoRoT- 
I's spectrum is well-reproduced by a 2460 K black-body to first order, which is interpreted 
to be due to either a high altitude layer that strongly absorbs stellar radiance, or an isother- 
mal region in the planetary atmosphere and that CoRoT-2 spectrum exhibits line emission 
from CO at 4.5 microns which could be caused by tidal- induced mass loss. The strength of 
the multi-wavelength approach is that it allows also to measure or place meaningful upper 
limits on the geometric albedo of a small set of transiting extrasolar planets. For example, 
the occultations of HAT-P-7b measured with the EPOXI, Kepler and Spitzer spacecrafts 
generated a broadband spectrum covering 0.35 — 8 /im. They ha ve also been used to iden - 
tify a set of atmospheric models that reproduce the observations ( IChristiansen et al.ll2010l ). 
The interpretation of this spectrum shows that it is consistent with an i nefficient day-night 
redistribution as seen from the Kepler phase curve ( iBorucki et al.ll2009l : ISpiegel &: Burrows 
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2010bl ). However, t he Kepler phase cu rve of this planet reveals a high night-side temperature 
at around 2600 K (jWelsh et al.ll2010l ). which appears to be at odds with knowing inefficient 
redistribution scenario a priori. 

In this paper, we combine occultation measurements obtained in the optical with Kepler 
and in the infra red with Warm-S yitzer to learn a bout the atmospheric properties of two 
HJs, Kepler-5b (IKoch et al.l I2OIOI ) and Kepler-6b (jPunham et al.ll2010l ). discovered using 
Kepler spacecraft. A search for occultations events was initially do ne using the two ffist 
Kepler quarters secured in long-cadence mode (IKipping fc Bakosll2010l ). In the present study, 
the photometry of the host stars was mearly continuously monitored by the Kepler space 
telescope allowing us to gather more than a year of observations and combine multiple 
occultations to improve the signal-to-noise. Thes e two planets have m asses and radii which 
are common amongst the known transiting HJs (ILatham et al.ll2010[ ). as shown in Tabled] 
Therefore, they may be representative of their class which make them good candidates to 
study, especially if we want to generalize our conclusions for comparative exoplanetology. 
We also obtained measurements of these targets in the near infrared with the Warm-Spitzer 
space telescope during occultations. Because of the high temperature of these HJs, our 
optical measurements are mainly due to thermal radiation rather than purely reflected light. 

We ffist describe the observations, time series and analysis for Warm-Spitzer data in 
Sect. 2 and Kepler data in Sect. 3 and finally discuss our results in Sect. 4. 



2. Warm-Spitzer observations and photometry 

This section describes the Warm-Spitzer observations and provides details of our ex- 
traction of the time series spanning each individual occultation event. 



2.1. Observations 



As Spitzer exhausted its cryogen of liquid coolant on 15 May 2009, only the ffist two 
channels, at 3.6 (c hannel 1) and 4.5 /im (channel 2), of the Infrared Array Camera (IRAC; 
Fazio et al.l ( l2004j )) are available in the post-cryogenic mission. Kepler-5 and Kepler-6 were 
both visited twice at each available IRAC bandpass, leading to a total 8 eclipses being se- 
cured as part of the 800 hours allocated to the program PID 60028 (PI: D. Charbonneau) . 
We present here the first observations from this program which contribute to characteriz- 
ing and vetting Kepler candidates. Each visit is secured in full frame mode (256 x 256) 
with an exposure time of 10.4 seconds, at 12 second cadence, leading to 2700 frames per 
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Kepler-5 (K=11.77) 10 hours observations and 2150 frames per Kepler-6 (K=11.71) 8 hours 
observations. The full set of Warm-Spitzer observations is presented in Table [2l 



2.2. Photometry 

We use the Basic Calibrated Data (BCD) frames produced by the standard IRAC cal- 
ibration pipeline for the photometric extraction. These files are corrected for dark current, 
flat-fielding, and detector non-linearity and converted into flux units. 

The first step of the photometric extraction consists of determining the centroid position 
of the stellar point spread function (PSF) using DAOPHOT-type Photometry Procedures, 
CNTRD, from the IDL Astronomy Library!^. We use the APER routine to perform an aperture 
photometry with a circular aperture of variable radius, using radii of 1.5 to 6 pixels, in 0.5 
steps. Finally, the propagated uncertainties are derived as a function of the aperture radius, 
and we adopt the aperture which provides the smallest errors. We find that the eclipse 
depths and errors vary only weakly with the aperture radius for all the light-curves analyzed 
in this project. In the case of Kepler-5, the optimal apertures are around 5 and 2.5 pixels 



er-6 has a com 



at 3.6 and 4.5 fim respectively. Kep 
and which is 3.8 magnitudes fainter (IDunham et al 



janio n at an angular separation of 4.1" 



20101). This star is located at 3.5 IRAC 



pixels from our target of interest. Therefore, we fix our circular aperture at a radius of 2.5 
pixels for all the visits of this star in order to minimize the contribution of the companion. 

The background level for each frame is determined by two methods. We use APER to 
measure the median value of the pixels inside an annulus centered on the star with inner 
and outer radii of 9 and 16 pixels, respectively. We also estimate the background by fitting a 
Gaussian to the central region of the histogram of counts from the full array. The center of 
the Gaussian fit is adopted as the residual background intensity. These two methods provide 
similar values; therefore, we use the background vlaues determined from an annulus. 

The contribution of the background to the total flux from the stars Kepler-5 and Kepler-6 
is low in both IRAC band-passes, from 0.2% to 0.55% depending of the photometric aperture 
size. Therefore, photometric errors are not dominated by fluctuations in the background. We 
note that the background is negative for all channel 1 observations and for one observation 
secured in channel 2 . The 10.4 seconds exposures yield a typical signal-to- noise ratio (S/N) 
of 200 and 150 per individual observation at 3.6 and 4.5 /xm respectively. 

After producing the photometric time series, we use a sliding median filter that compares 



^http : //idlastro . gsf c . nasa.gov/homepage . html 
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the 20 preceding and 20 following photometric measurements and centroid positioning to 
identify and trim outliers greater than 5 a. This process removes measurements affected 
by transient hot pixels and inaccurate centroid determination. In this way, we discarded 
approximately 2% photometric points from all the observations. We also discarded the first 
half hour of observations, corresponding to around a hundred frames. These frames exhibit 
an anomalously large pointing drift, most likely due to settling of the telescope at the new 
position. The final number of photometric measurements used for each observation is given 
in Table H 



2.3. Occultation amplitudes and associated errors 

To measure the occultation depths and their uncertainties we model the light curves 
with 4 parameters: the occultation depth, d, the orbital semi-major axis to stellar radius 
ratio (system scale), a/i?^, the impact parameter, 6, and the time of mid transit, Tc. We use 
the IDL transit routine OCCULTSMALL, developed by Mandel & Agol (2002), to model the 
light curve. Kepler-5b and Kepler-6b have both well defined transit and stellar parameters 
constrained by the transits from the Kepler light-curves (Koch et al. 2010, Dunham et 
al. 2010). Thus, we adopt these parameters, in particular the inclination and the scale of 
the system, by fixing a/i?^ and h to their nominal values in our eclipsing model. We also 
fixed the mid-eclipse times to the values we derived from the Kepler lightcurves (see Sect. [3] 
below). Only the depth of the occultations are allowed to vary when fitting the Warm-Spitzer 
observations. 

The Warm-Spitzer /\R AC photometry is known to be systematically affected by the 
so-called pixel-phase effect, which is due to the combination of pointing jitter and intra- 
pixel sensitivity (see, e.g., Charbonneau et al. 2005, 2008, Reach et al. 2006). This effect 
corresponds to oscillations in the measured raw light curve with an approximate period of 
70 minutes and an amplitude of 2% peak-to-peak. This artefact has to be corrected to 
properly extract the eclipse depths and errors. To correct the raw light curve for this intra- 
pixel sensitivity effect, we use the centroid position of the target on the detector and its 
variations as function of time. We use a quadratic function of the position and a linear 
function of time F^orr = F[Ki{x -xq) + K2{x -xq)^ + Ks{y -yo) + K^iy -yof + K5 + K^xt], 
where F and F^orr are the fluxes of the star before and after the pixel-phase effect correction, 
and {x — Xq) and {y — yo) define the position in pixel of the source centroid on the detector 
with respect to the pixel pointing position, located at [xc^/o], and the constants Ki are the 
free parameters to adjust. 
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We use the MPFIT packag^ to perform a Levenberg-Marquardt least-squares fit of the 
transit model to observations. The best-fit model is computed over the whole parameter 
space [d, Ki). The baseline function described above is combined with the transit light 
curve function so that the fit is constrained by 8 parameters (2 for the transit model, 2 for 
the linear baseline, and 4 for the pixel phase effect). 

We use four methods to determine the centroid position of a stellar PSF and test the 
its impact on the final results. In the first method, we fit Gaussians to the marginal x and 
y sums using GCNTRD. The method we apply consists of computing the centroid of the stars 
using a derivative search with CNTRD also part of the standard IDL Astronomy Library. As 
for the third method, we calculate the center-of-light of the star within an circular aperture 
with a radius of 3.0 pixels to approximate the center of the star. Finally, the fourth method 
uses a symmetric 2D Gaussian fit with a fixed wid t h to a 3x3 pixel sub-array to approximate 



the center of the star, as suggested by lAgol et al.l (120101 ). 



All these methods allow us to measure the centroid posi tion of a stellar i mage with a 



different level of accuracy and precision as previously noted by lAgol et al.l (120101 ). For all the 
methods tested, we find that the centroid position varies by less than 20% of a pixel during 
a complete observation and that can be determined to a precision of a hundredth of a pixel. 
We measure the eclipse depths and associated errors for all the four methods following the 
algorithm we describe below. Although each method provides slightly different results for 
the centroid determination, it yields an eclipse depth consistent to within 1 a and similar 
errors. In this paper we use the CNTRD program as it produces the smallest reduced . 

As shown by Pont et al. (2006), the existence of low- frequency correlated noise (red 
noise) between different exposures must be considered to obtain a realistic estimation of the 
uncertainties. To obtain an estimate of the systematic errors in our observations we use the 
permutation of the residuals method, known as "prayer-bead" (Moutou et al. 2004, Gillon 
et al. 2006, 2007) and derived the covariance from the residuals of the light curve. In this 
method, the residuals of the initial fit are shifted systematically and sequentially by one 
frame, and then added to the transit light curve model before fitting again. The error on 
each photometric point is the same and is set to the rms of the residuals of the first best-fit 
obtained. We find that the final rms values are 5 to 20% larger than the predicted photon 
noise depending on the target and the bandpass. 

We produce as many shifts and fits of transit light curves as the number of photometric 
measurements to determine the statistical and systematical errors for the adjusted parame- 
ters. We set our uncertainties equal to the range of values containing 68% of the points in 



^http : //cow . physics . wise . edu/^craigm/idl/idl . html 
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the distribution in a symmetric range about the median for a given parameter. We check 
that these values are close to the standard deviation of each nearly Gaussian distributions. 
We fit for the eclipse depth fixing the mid-eclipses time to the values we derived from the 
Kepler observations as describe in the following section. The best-fits of the occultation 
depths their error bars are listed in Table |2l 



3. Analysis of the Kepler light curves 



This section describes the Kepler observations gathered almost continuously during 
nearly 14 months and provides details on the analysis of the different occultations. 

The Kepler pass -band spans 437 to 897 nm, with a central wavelength rou ghly equiv- 



alent to the R-band ( iKoch et al.ll2010t iBatalha et al.ll2010l : iBryson et al. 



Kepler-5 and Kepler -6 have Kepler magnitudes (Kp) of 13.37 and 13.30 ( iKoch et al.l 12010 



2010). The stars 



Dunham et al.l 120101 ). The light of the companion star close-by Kepler-6 as discussed above 
in Sect. [2] is included in the Kepler aperture. This has the effect o f dilut ing the depth of 
the eclipse signal by only a few percent as shown by iDunham et al.l (|2010[ ) , well below the 
precision we obtain (see below). 

We use the Kepler science data of Kepler-5 and Kepler-6 from Quarter to 5 (Q0-Q5). 



Thes e observations have been reduced and detrended by the Kepler pipeline ( 1 Jenkins et al. 



20101). The pipehne produces both calibrated light curves (PA data) for individual analysis 
and corrected light curves (PDC) which are used to search for transits. We used both datasets 
and checked that the conclusions remain the same. This paper presents results which are 
measured from P DC data. They consi s t of long cadence int egration time (30 minutes) for 
Quarters and 1 ICaldwell et al.l (l2010l) : iJenkins et al.l ( 120101 ) and short cadence (1 minute) 
for Quarters 2 to 4 lGilliland et al.l ( 120101 ). We acquired the data in a uniform manner at short 
cadence between BJD 2454998 and 2455371 (UT 2009 June 15 - 2010 June 23) during 373 
days of observations (Q2-Q5). The pipeline provides time series with times in barycentric 
corrected Julian days, and flux in electrons. 

The quality of the photometric time series around each expected eclipse event is checked 
visually in a first step. This step allows to flag and remove eclipses that we consider not 
satisfactory for various reasons (artefacts, Spacecraft-related events occasionally result in loss 
and subsequent reacquisition of fine guidance, etc.). We reject 20 of the 103 occultation 
events observed for Kepler-5, and 12 of these events being amongst 113 for Kepler-6. The 
numbers of occultations that we keep per target and per observing mode are presented in 
Table 131 We then identify and trim outliers greater than 5 a using a sliding median filter that 
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compares the 20 preceding and 20 following photometric measurements. This process rejects 
a very small number of photometric measurements (less than 0.01%) compared to the total 
measurement number (102,300 for Kepler-5 and 114,000 for Kepler-6). At this point we check 
that the rms scatter is constant over the whole observational period. We fit each individual 
eclipses between the planetary orbital phase 0.38 to 0.62 using the occultation model (see 
description in Sect. 12. 3p assuming constant ephemerides and null limb darkening coefficients. 
All the transit parameters but the occultation depths and the mid-occultation times are 
set to fixe values. A linear function of time represents the baseline. Each individual fitted 
light-curves is then normalized to its local baseline in order to produce a set of normalized 
eclipse light curves. We finally phase fold and combine all the normalize light curves to 
produce four photometric time series, one for each observing mode (cadence) and target. 
We present the normalized, folded and combined light curves obtained at short cadence for 
both targets in Figure E] from which we measure the occultation depths as described in the 
following paragraphs. 

We first search for the occultation events in the combined and normalized occultation 
light curves. This is done by evaluating the depths of Levenberg-Marquardt fits to models 
with the same overall shape as the planetary transits (same impact parameters and system 
scales), but at different phases of the orbital period. The best fits as function of the orbital 
phases are plotted in Fig. HI and the maximum depths and best fits are found very close to 
the orbital phase of 0.5 as expected for a circular orbit. This suggests that we detect the 
occultations for both targets in the Kepler observations. 

We now estimate the significance of these detections by measuring the occultations 
depths, ephemerides and associated errors for each observing cadence mode and targets. We 
estimate the parameter values and errors using a bootstrap Monte-Carlo analysis. In this 
method, we first find the best-fit to the occultation curves and produce a set of photometric 
residuals. We compute the rms of the residuals and create a new simulated data set by adding 
the errors, chosen randomly from a Gaussian distribution having the rms for its FWHM, to 
the best-fit. We fit the new simulated data set and repeat this process for 10^ trials. We 
compute the median value and the standard deviation of the Gaussian distributions to derive 
the occultation depth and the central phase errors. The depths and mid-occultations as well 
as measurement significances are plotted as function of the orbital phases in Figure M and 
presented in Table [31 We weighted average the best-fitted depths and times to obtain one 
single value per target. We find that both planet have similar occultation depths of 21 ± 6 
and 22 ± 7 ppm for Kepler-5b and Kepler-6b respectively. 
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4. Results 

We have shown that decreases in the fluxes of the Kepler and Warm-Spitzer light-curves 
correlated with the expected phases of occultations are detected in all band-passes at 3 a 
level. We thus interpreted these signals as planetary occultations detected in the optical and 
in the near infrared wavelengths. 



4.1. Constraining the orbital eccentricity 

The current upper limit on the orbital eccentricity e from radial velocity measurements 
is consistent with zero for both targets (Koch et al. 2010, Dunham et al. 2010). An 
i ndependent analysis of the two first Kepler quarters (QO and Ql) have already been applied 



( iKipping fc Bakosll2010[ ). It uses the transits, eclipses and radial velocity measurements to 
confirm that no eccentricity is found for Kepler-5 and to report a marginal eccentricity at 
2 a level for Kepler-6. 

We measure the mid-occultation timing offset from both Kepler and Warm-Spitzer 
observations. The determination of the timing of the secondary eclipse constrains the planet's 
orbital eccentricity. A non-zero value of e could produce a measurable shift in mid-eclipses 
times. The Warm-Spitzer observations are affected by systematics (correlated noise) which 
prevent us from deriving timing measurements with a better precision than the precision 
we derive from Kepler measurements. Indeed, we measure the mid-occultation timing offset 
from Kepler with uncertainties of nearly 8 minutes for both targets. Our estimate for 
the best-fit timing offset translates to a constraint on e and the argument of pericenter u. 
The timing is used to constrain the ecos(a;) at 3 cr level. For Kepler-5b we find ecos(u;) = 
-0.025 ±0.005, and |ecos(w)| < 0.06 to 3 a. For Kepler-6b we find ecos(u;) = -0.01 ±0.005, 
and |ecos((X')| < 0.035 to 3 a. These upper limits imply that the orbits of these objects are 
nearly circular unless the line of sight is aligned with the planet's major axes, i.e. the 
argument of periapse u is close to 90° or 270°. Assuming that both planets are on a circular 
orbit and accounting for the time taken by the light to cross the entire orbits, we would expect 
that the mid-secondary eclipses are delayed by 50 seconds and 45 seconds for Kepler-5 and 
Kepler-6 respectively, which translate into orbital phase of 0.50016 for both targets. 



4.2. Atmospheric considerations 



The occultation depths measured in each bandpass are combined and are turned into 
an emergent spectrum for each planet. The observed flux of the planet in each bandpass 



- 11 - 



corresponds to the sum of the reflected hght, the thermal emission of the incident stellar flux, 
and the interior flux from the planet itself (e.g., internal heat or emission due to tidal forces). 
As a first step, we assume that the thermal emission from the planet itself is negligible. We 
use the chromatic information to distinguish among the different components of the observed 
planetary flux. 

If we were to explain the occultation detections in the Kepler bandpass as reflected light 
for both objects, it would imply that we are indeed measuring their geometric albedos. One 

can estimate the reflected light by the planet as -p- = (j^^ where Ag is the geometric 
albedo in the Kepler bandpass. With this assumption in mind, and with the occultation 
depths of 21 ±6 ppm and 22 ±7 ppm for Kepler-5b and Kepler-6b respectively (see Table |3]), 
we find that both planets have low geometric albedos in the Kepler band-pass. Thus, we 
obtain Ag = 0.12 ± 0.04 and Ag = 0.11 ± 0.04 for Kepler-5b and Kepler-6b respectively. 
Interestingly, we can infer the Bond- albedo Ab from Ag. The Bond albedo is the fraction 
of the bolometric incident radiation that is scattered back out into space at all phase angle. 
Since we have no phase curves information for both objects, we assume a Lambertian criteria 
where A^ < 1.5 x Ag, we which leads to A^ < 0.17 at 1 cr level for both object. 

The equilibrium temperature is derived from 

Te, = T,(i?»l/2[/(l-AB)]l/4 (1) 

which depends on the Bond albedo Ab and the re-distribution factor / which accounts for the 
efficiency of the transport of energy from the day to the night side of the planet. / can vary 
between 1/4 for an extremely efficient redistribution (isothermal emission at every location 
of the planet) and higher values for an inefficient redistribution, implying big differences 
in the temperatures between the day and the night sides of the planet. Nevertheless, the 
detection of the occultations in the Kepler bandpass may indicate that we measure the 
thermal emission of these planets at these wavelengths, as expected for the temperature 
regime of HJs. Therefore, the Bond-albedo could be indeed well below 0.17. Notably, as 
Kepler spacecraft continues to monitor the photometry of Kepler-5b and Kepler-6b, the 
observations will provide better constraints on occultations and may reveal the phase curves 
which would allow to fully constrain the Bond-albedo of these planets. 

We estimate the thermal component of the planet's emission in the two Warm-Spitzer 
band-passes from the occultation depths measured at 3.6 and 4.5 /xm (see Table To do 
so, we assume that the planetary emission is well reproduced by a black-body spectrum and 
translate the measured depth of the secondary eclipse into brightness temperatures. We 



use the PHOENIX atmospheric code (iHauschildt et al.l 119991 ) to produce theoretical stellar 



models for the star, with stellar temperature of Teg = 6297 ± 60 for Kepler-5b (Koch et 



- 12 - 



al. (2010) Teff = 5647 ± 44 for Kepler-6b (Dunham et al. 2010). Taking the Warm-Spitzer 
spectral response function into account, the ratio of areas of the star and the planet and 
the stellar spectra, we derive the brightness temperatures that best fit the observed eclipse 
depths measured in the two IRAC bandpasses. The brightness temperature calculated this 
way resulted in Tspitzer=1930±100 K and Tspitzer=1660± 120 K for Kepler-5b and Kepler-6b 
respectively. If we further assume that the planet is in thermal equilibrium and has a zero 
Bond albedo, these temperatures favor low values of the re-distribution factor /, i.e. the 
planet day-side re-radiate efficiently the incoming stellar energy flux. 

The brightness temperatures we derive using the Warm-Spitzer band-passes are larger 
than the equilibrium temperatures for both targets (see Table. [1]). This may indicates that 
our assumption that these planets behave as black-bodies is most probably wrong, simi- 
lar to the giant planets of our solar system. Furthermore, a more robust determination 
of requires deta iled model computation because of the wide range of parameter space 



( ISeager et al.l 120051 ). We discuss below our results from two different types of atmospheric 



models. 



In the first approach, we use the atmospheric retrieval method of lMadhusudhan fc Seager 



( 120091 ) to derive the temperature structure and composition of each system, given the data. 
The model involves 1-D line-by-line radiative transfer, with parametric temperature struc- 
ture and composition, and includes the major molecular and continuum opacity sources, 
along with constraints of LTE, hydrostatic equilibrium, and global energy balance. This 
modeling approach allows one to compute large ensembles of models (~ 10^), and explore 
the parameter space of molecular compositions and temperature structure in search of the 
best-fitting models. In the present work, the dominant sources of opacity included are: H2O, 
CO, CH4, CO2, NH3, and H2-H2 CIA absorptions, in the infrared, and TiO, VO, Na, K, and 
Rayleigh scatteri ng, in the visible. In the prese nt context, the number of model parameters 



is = 10 — 15 (IMadhusudhan fc Seagerll2009l ). Thus, the limited number of available ob 



servations (Aobs = 3) imply a substantial degeneracy in solutions. Consequently, our goal 
with the present data is to find a nominal set of solutions, as opposed to finding unique fits. 

Models fitting the observations of Kepler-5b and Kepler-6b are shown in Figure |Hi We 
find that observations of Kepler-5b can be explained to good precision by a wide range in com- 
position, including models with equihbrium chemistry assuming solar abundances. However, 
current observations tentatively rule out a thermal inversion in this syste m. At the temper 



ature s of Kepler-5b the atmosphere is expected to be abundant in CO ( jBurrows fc Sharp 



19991 ). which has a strong feature in the 4.5 /xm bandpass. A thermal inversion would there- 
fore naturally predict a high fiux in this channel over the 3.6 /xm channel, as opposed to 
the observed fiuxes which are similar between the two channels. Instead, the CO absorp- 
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tion feature caused by the lack of a thermal inversion explains the observations very well. 
Additionally, the Warm-Spitzer and Kepler data together require that the incident energy 
is mostly re-radiated on the dayside, with a low albedo and a small fraction (/) of dayside 
energy redistributed to the night side. Assuming zero albedo, the model shown has / = 0.12. 

The observations of Kepler-6b allow for the possibility of a thermal inversion in its 
atmosphere. As shown in Figure [H the data can be explained to within the l-cr errors 
by models with and without thermal inversions. As explained above, the higher flux in 
the 4.5 /im channel over the 3.6 /xm channel can be interpreted as a sign of an inversion. 
However, given the degrees of freedom allowed by the molecular abundances, the limited 
data can be fit just as well without inversions (see e.g. Madhusudhan & Seager, 2010b). A 
wide range of abundances can fit the data, including those close (within a factor of ten) to 
equilibrium chemistry assuming solar abundances. However, as with Kepler-5b the Warm- 
Spitzer and Kepler data together place stringent constraints on the energy budget on the 
dayside, requiring low albedos and low redistribution to the night side. 

The second mod el we consider to compare our data to is the HJ atmospheric model of 
Fortney et al.l (120081 ). Our motivation here is that, given our limited wavelength coverage, 
it is important to compare our best-fit planetary-to-star flux ratio values with physically 
motivated models. We aim at broadly distinguishing these atmospheres between different 
classes of models, given that there are no prior constraints on the basic composition and 
structure. 



The atmospheric spectrum calculations are performed for ID atmospheric pressure- 



Lodders 


( 


2002 


); 



of radiative transfer and chemical e quilibrium of neutral species. The opacity database is 



described in iFreedman et al.l (120081 ). So far, this model has been used to generate P-T 
proflles for a variety of close-in planets iFortney et al.l (l2005l . |2006| . 120081 ). The models are 
calculated for various dayside-to-nightside energy redistribution parameter (/) and allow for 
the presence of TiO at high altitude, playing the role of absorber which likely lead to an 
inversion of the T-P proflle. 

We compare the data to model predictions and select models with the best reduced . 
We note that this is not a flt involving adjustable parameters. All models for both targets 
show that the dayside is re-radiating the stellar flux efficiently, favoring low values for /. 
Since our observations are weakly constraining, model comparisons for Kepler-6b (Figure [H]) 
suggest that both inverted and non-inverted atmospheric temperature proflles can reproduce 
the data. However this is not the case for Kepler-5b where the models show no inversions. 
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For HJs such as Kepler-5 and Kepler-6, both scattered stellar light and planetary 
thermal emission could contri bute to the planet emergent flux in the Kepler bandpass 



(iLopez- Morales fc Seagerl 120071 ) . The ratio of scattered starlight to thermal emission de- 
pends on the atmospheric composition. These objects could be too hot for condensates. 
Without a reflective condensate layer, photons in the Kepler bandpass are absorbed before 
being scattered. Therefore, the albedo of Kepler-5 and Kepler-6 are likely to be low, as 
observed here. 



Knutson et al.l (120101 ) shows that there could be a correlation between the host star 
activity level and the thermal inversion of the planetary atmosphere. Effectively, the strong 
XEUV irradiation from the active stellar host of HJs could deplete the atmosphere of chem- 
ical species responsible for producing inversions. The Ca II H & K line strengths are a good 
indicator of the stellar activity. We obtain Keck HIRES spectra for Kepler-5 and Kepler-6, 
and estimate the line strengths of S'hk=0.14 and logi?HK="4.97 for Kepler-5 and S'hk=0.16 
and log i?HK="4.98 for Kepler-6. Both star are moderately quiet sub-giants , with log-R^y 



falling in-between cases for thermal inversion or no inversio n (see discussion in iKnutson et al. 



(120 lOl )). We also evaluate the empirical index defined by iKnutson et al.l (120101) which pro 



vides a way to distinguish between the different hot- Jupiter atmospheres. IKnutson et al. 



(I2OIOI ) proposes an index value that could be correlated with the presence on a thermal 
inversion. Using the same definition, we find that the index = —0.023 ± 0.023 and index = 
0.040 ± 0.033 for Kepler-5 and Kepler-6 respectively. These values suggest that the atmo- 
sphere of Kepler-6b is consistent with a weak thermal inversion and the one of Kepler-5b 
with a non-inverted profile, which are both in agreement with the results of our present 
study. 

As Kepler spacecraft continues to monitor the photometry of Kepler-5b and Kepler-6b, 
the future quarters of observations will be used to improve the signal-to-noise. This will 
provide better constraints on the occultations and may reveal the phase curves which would 
allow to derive a more accurate Bond albedo and thus provide better constrain the energy 
budget of these planetary atmospheres. 



We would like to thank Jacob Bean and Brice-Olivier Demory for a variety of useful 
discussions. This work is based on observations made with the Spitzer Space Telescope, 
which is operated by the Jet Propulsion Laboratory, California Institute of Technology un- 
der a contract with NASA. Support for this work was provided by NASA through an award 
issued by JPL/Caltech. This work is also based on observations made with Kepler which 
was competitively selected as the tenth Discovery mission. Funding for this mission is pro- 
vided by NASA's Science Mission Directorate. The authors would like to thank the many 
people who generously gave so much their time to make this Mission a success. Some of 
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Table 1. System Parameters for Kepler-5b for Kepler-6b obtained from the discovery 

papers. 



Parameter 



Value 



Kepler-Sb from Koch et al.(2010) 



Orbital period P (d) 
Mid-transit time E (HJD) 
Scaled semi-major axis a/R^, 
Scaled planet radius i?p/R« 
Impact parameter b = a cos i/R* 
Orbital inclination i (dog) 
Orbital eccentricity e 

Stellar parameters 
Effective temperature T^g (K) 
Mass M,(M0) 
Radius R*{Rq) 

Planetary parameters 
Mass Mp (Mj) 
Radius Rp (Rj, equatorial) 
Density pp (gcm^''') 
Orbital semi-major axis a (AU) 
Equilibrium temperature Teq (K) 



3.548460 ±0.000032 
2454955.90122 ± 0.00021 
6.06 ±0.14 
n8195+° °°°3° 

U.UOi90_Q QQ047 

Qqq+0.061 

86?3±0.5 
< 0.024 



6297 ± 60 

1 Q7/1+0.040 

1 7QQ+0.043 
1. (»o_g Qg2 



2.114 

1.431 



0.056 
-0.059 
0.041 



-0.052 

0.894 ± 0.079 
0.05064 ± 0.00070 
1868 ± 284 



Kepler-6b from Dunham et al.(2010) 



Orbital period P (d) 
Mid-transit time E (HJD) 
Scaled semi-major axis a/R^ 
Scaled planet radius Rp/R* 
Impact parameter b = a cos i/Ri, 
Orbital inclination i 
Orbital eccentricity e 

Stellar parameters 
Effective temperature T^g (K) 
Mass M*(M0) 
Radius R*{Rq) 

Planetary parameters 
Mass Mp (Mj) 
Radius Rp {Rj, equatorial) 
Density pp (g cm~^) 
Orbital semi-major axis a (AU) 
Equilibrium temperature Teq (K) 



3.234723 ± 0.000017 
2454954.48636 ± 0.00014 

7 05+"-" 

0.09829iHSS60 
0.398l°:°f, 
86?8±0.3 
(adopted) 



5647 ± 44 



1.209 



-0.044 



-0.038 
1 391+°°l'' 



0.669 
1.323 



0.025 
0.030 

+0.026 
0.029 
0.018 



0.352_(, „22 
n n45fi7+"''-"-"-'^^ 

U.U<iOD/_Q QgQ4g 

1500 ± 200 
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Table 2. Warm-Spitzer ohseivations. 



Target Visit Wavelength Obs. Date (UT) Select, points 



Depth 



Weighted. Avg. depth Bright. T (K) 



Kepler- 5 
Kepler- 5 
Kepler-5 
Kepler-5 



3.6 
3.6 
4.5 
4.5 



2009- 12-26 

2010- 06-16 
2010-01-05 
2010-07-12 



2455 
2465 
2540 
2242 



0.091 ± 0.022% 
0.120 ± 0.026% 
0.102 ±0.023% 
0.111 ±0.021% 



0.103 ±0.017% 
0.107 ±0.015% 



1900 ±110 



1770 ± 100 



Kepler-6 
Kepler-6 
Kepler-6 
Kepler-6 



3.6 
3.6 
4.5 
4.5 



2009- 12-30 

2010- 06-16 
2010-01-02 
2010-01-13 



2013 
2012 
2022 
1993 



0.108 ± 0.048% 
0.051 ± 0.032% 
0.123 ±0.027% 
0.180 ±0.027% 



0.069 ± 0.027% 
0.151 ±0.019% 



1320 ± 250 
1700 ± 120 



Table 3. Kepler observations. 



Target 


Quarters 


Nbr of eclipses 


Observing mode 


Fp/F^, (ppm) 


Weighted. Avg. Fp/Ft (ppm) 




Kcplcr-5 
Kepler-5 


to 1 
2 to 4 


11 
83 


Long cadence 
Short cadence 


25 ± 19 
2ll| 


21 ±6 


0.12 ±0.04 


Kepler-6 
Kepler-6 


to 1 
2 to 4 


12 
101 


Long cadence 
Short cadence 


27± 19 

21+* 


22 ±7 


0.11 ±0.04 
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Kepler-5b Spltzer/IRAC Kepler-6b Spltzer/IRAC 




-6 -4 -2 2 4 -6 -4 -2 2 4 

Time [Hours] Time [Hours] 



Fig. 1. — Warm-Spitzer raw occultation light-curves observed in the two currently available 
IRAC band-passes at 3.6 and 4.5 /im, showing the two visits per bandpass (Table. E]). Plot 
for Kepler-5b is on the left side and Kepler-6b is on the right side. The data are binned in 
12 minute intervals and offset in flux for display. The red solid line corresponds to the best 
fit models which include the time and position decorrelation as well as the model for the 
planetary occultation (see details in Sect. 12. 3p . 



Kepler-5b Spltzer/IRAC Kepler-6b Spltzer/IRAC 




Time [Hours] Time [Hours] 



Fig. 2. — De-correlated, normalized and combined per bandpass Warm-Spitzer occultation 
light-curves. The best fit are over-plotted in red solid line. The best fit depths with their 
associated errors are provided in Table. |2] 

the data presented herein were obtained at the W.M. Keck Observatory, which is operated 
as a scientific partnership among the California Institute of Technology, the University of 
California and the National Aeronautics and Space Administration. The Observatory was 
made possible by the generous financial support of the W.M. Keck Foundation. 
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Kepler 



1.00002 



1.00000 



0.99998 



0.99996 




Kepler-6b 



0.35 0.40 



0.45 0.50 0.55 
Orbital phase 



0.60 0.65 



1.00002 



- 1.00000 



0.99998 



0.99996 



0.35 0.40 



0.45 0.50 0.55 0.60 
Orbital phase 



Fig. 3. — Phase folded normalized Kepler occultation light-curves. Plot for Kepler-5b is 
on the left side and Kepler-6b is on the right side. The observation are the black point 
with their one-sigma error-bars. The data are binned by approximately 40 minutes (0.01 in 
phase). The best fit are over-plotted with the red solid line. 
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Fig. 4. — Delta from the best fit of the depth of a model as function of the mid-occultation 
time from zero (phase 0.5). Plot for Kepler-5b is on the left side and Kepler-6b is on the right 
side. The shape and duration of the occultation model are set by the transiting parameters. 
The vertical dashed lines show the transit durations and around zeros which are highlighted 
with vertical dotted lines. Mid-occultations are expected at phase 0.5 for non eccentric 
orbits. The horizontal dotted lines indicates the limits of 1, 2 and 3 a from the best fit. The 
maximum occultation depth is centered near phase 0.5 which confirm the detection of the 
occultations event. 
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Fig. 7. — The Kepler and Warm-Spitzer occulta tions measured for Kepler-5b and Kepler- 
6b compared to hot- Jupiter atmospheric models (IMadhusudhan &: Seagerll2009l ). The black 
points associated with their one-sigma error bars are the measurements. The solid lines 
correspond to the best fit models. The green and red points (without error bars) are the 
bandpass-integrated model values. The dashed line show black-body ratios for Warm-Spitzer 
measurements at various temperatures (see details in Sect. 14.21) . The insets are magnifications 
of the optical region of the spectra to reveal the Kepler measurements. 
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Fig. 8. — The temperature-pressure profiles for Kepler-5b and Kepler-6b corresponding to 
tlie model presented in Figure |H1 The profile with a thermal inversion is over-plotted in red. 
Both profiles fits well Kepler-6b observations. 
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Fig. 9. — Day s ide planet-to-star fl ux ratios as function of the wavelength for three atmo- 
spheric models (IFortney et al.ll2008l ) per target. The orange diamonds and their error bars 
correspond to the Kepler and Warm-Spitzer observations. The blue models represent non- 
inverted atmospheres (no TiO) and no redistribution of the energy to the planetary night 
side. The green models represent inverted atmospheres (with TiO) and no redistribution of 
the energy to the planetary night side. The red models represent inverted atmospheres (with 
TiO) and with full redistribution of the energy to the planetary night side. Both targets 
are best fitted with non-inverted models which include dayside redistribution, with also an 
equivalent good fit for an inverted atmosphere for Kepler-6b. 



